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Abstract. We give a simple proof of LP boundedness of iterated commutators of Riesz 
transforms and a product BMO function. We use a representation of the Riesz transforms 
by means of simple dyadic operators - dyadic shifts - which in turn reduces the estimate 
quickly to paraproduct estimates. 



I. Introduction 

It was shown by the authors in [14] that product BMO of S.-Y. A Chang and R. Fefferman, 
defined on the space M. dl <g) • • ■ <g> M. dt , can be characterized by the iterated commutators of 
Riesz transforms in the multi-parameter setting. This extended a classical one-parameter 
result of R. Coifman, R. Rochberg, and G. Weiss [3], and at the same time extended the 
work of M. Lacey and S. Ferguson [7] and M. Lacey and E. Terwilleger [13], which treats the 
case of iterated commutators with Hilbert transforms. 

In the multi parameter setting, one seeks to prove an inequality of the type 

(1) (2) 

Memo < \\[T 1 ,[T 2 ,...[T n ,M b ]...]]\\L^ < Memo 

where Tj are either Riesz or Hilbert transforms acting in the ith variable and M& is multipli- 
cation by b. We refer to (1) as a lower bound and (2) as an upper bound. 

In the case of the Hilbert transforms, the upper bound is not hard to prove, by appealing to 
analyticity even in the multi parameter setting, see [8] , whereas the lower bound is extremely 
difficult and makes use of the corresponding upper bound. Similar considerations for the 
Riesz transforms fail because of a lack of analyticity. The upper bound is already difficult, 
even in the case of one parameter, and the proof of the lower bound is non-trivial. Coifman, 
Rochberg and Weiss [3] used a sharp function method to prove their one-parameter analog of 
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Theorem 1.1. This method, versatile as it is in the one-parameter case, has no clear analog in 
the product setting; see Lacey [12] for a discussion of this in a similar context. The methods 
of [3] do not extend to the multi parameter setting. 

In [14] we find it necessary to prove an upper bound for more generalized smooth cone 
projection operators to simulate the analyticity we do not have. Such operators, even in 
the one parameter case, are not included in the work in [3]. The upper bound for the Riesz 
transforms in particular allow for a relatively simple proof that extends naturally to the multi 
parameter setting. In this note, we present this proof. Here is the theorem: 

1.1. Theorem. 

\\[R ldl ,[R 2>j2 ,...[R n>jn ,M b ]...]]\\ L ^ L 2 < \\b\\ BMO 

Here R s j denotes the jth Riesz transform acting on the sth variable on M ds and Mf, denotes 
multiplication by the symbol b. The BMO norm used is Cheng Fefferman product BMO. 

The proof given here is rather different from that of Coifman, Rochberg and Weiss. The 
proof of the upper bound presented here follows from the decomposition of the commutators 
into a sum of simpler terms. The decomposition we allude to is naturally described in terms 
of 'Haar shifts,' a method of proof used by Petermichl [18] to address a subtle question 
concerning the Hilbert transform in a context with matrix valued weights and then extended 
by Petermichl, Treil and Volberg in [19] to the Riesz transforms. This method reduces the 
main issue to paraproduct estimates quickly. The paraproducts that arise are of multi- 
parameter form. The specific result needed is Theorem 5.1 below. This result is due to 
Journe [9]; more recent discussions of paraproducts are in [11,16,17]. 

2. The Upper Bound 

def 

Let M;, if = bip be the operator of pointwise multiplication by a function b. For Schwartz 
functions / on M d , let Rj / denote the jth Riesz transform of /, for 1 < j < d. 

We are concerned with product spaces M. d = M dl (g>- ■ ■®M a! ' for vectors d = (d\, . . . , d t ) G N*. 
For Schwartz functions b, f on M. d , and for a vector j = (ji, . . . ,jt) with 1 < j s < d s for 
s — 1, . . . , t we consider the family of commutators 

(2-1) Cj(b,f)(x) d = f [•■•[[M 6 ,R 1)il ],R 2 , j2 ] ) ...](/)(x) 

where K s j denotes the j'th Riesz transform acting on M. da . 

2.2. Theorem. We have the estimates below, valid for 1 < p < oo. 

(2-3) l|Cj(6,^)|| p < ||&||bmo|MIp- 

By BMO, we mean Chang-Fefferman BMO. 
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We start in the one-parameter setting. We will use dilation and translation operators on 
R d 

(3.1) Tr y f(x) = f(x-y), y G R d , 

(3.2) Dili p) f{x) = a- d/p f(x/a) , a>0,0<p<oo. 
These will also be applied to sets, in an obvious fashion, in the case of p = oo. 

By the canonical dyadic grid in M. d we mean the collection of cubes 

V cncl = {j2 k + [0, 2 k ) d : j e z d , k e Z} 

By a dyadic grid we mean any of the collections 

V d = Vf = {Dil[ oo) Ti y I : / G V cncl } , 1 < t < 2 , y e M d . 

Note that T> cnc \ = T> d '°. If we emphasize the role of dimension, we refer to these collections 
as d- dimensional dyadic grids. 

Haar functions on M d are now described. For e G {0, 1}, set 

lO def n n ,i def ^ 

h - - 1 hio) + 1 [oi)' h = Vi|)" 

Here, we put the superscript to denote that 'the function has mean 0,' while a superscript 
1 denotes that 'the function is an L 2 normalized indicator function.' In one dimension, for 
an interval /, set 



e 



Of course for any choice of one dimensional grid T>, the collections of functions {h® : / e T>} 
form a Haar basis for L P (M). 

Let Sig d = f {0, l} d — {1}, which we refer to as signatures. In d dimensions, for a cube Q 
with side J, i.e., Q = I x • • • x J, and a choice of e G Sig d , set 

d 

h £ Q (x lj ...,x d ) = IJ/ijfo)- 

J"=l 

It is then the case that the collection of functions 

Haar£> d = {h e Q : Q e V d , e e Sig d } 

form a Haar basis for L p (M d ) for any choice of ci-dimensional dyadic grid V d . Here, we are 
using the notation 1 = (1,...,1). While we exclude the superscript 1 here, they play a role 
in the theory of paraproducts. 
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We will use these bases in the tensor product setting. Thus, for a vector d = (di, . . . , d t ), 
and 1 < s < t, let X> ds be a choice of d s dimensional dyadic grid, and let 

Also, let Sig^* = f {e = (ex, ...,£*) : e s G Sig ds }. Note that each e s is a vector, and so e is 
a 'vector of vectors'. For a rectangle R — Qi x ■ ■ -Q t , i.e. a product of cubes of possibly 
different dimensions, and a choice of vectors e G Sig^- set 

t 

h R (x 1 ,...,x t ) =J[h^ s (x s ) 

s=l 

These are the appropriate functions and bases to analyze multi-parameter paraproducts and 
commutators. 

Let 

Haarx> d - = f {h £ R : ReVj, EE Sig d -} . 
This is a basis in L p (IR d ), where we use the notation 

R d d i R <ii 8 ... 8 R * 

to emphasize that we are in a tensor product setting. 

4. Chang-Fefferman BMO 

We describe the elements of product Hardy space theory, as developed by S.-Y. Chang 
and R. Fefferman [1,2,4-6] as well as Journe [9,10]. By this, we mean the Hardy spaces 
associated with domains like ®^ =1 M ds . 

4.1. Remark. The (real) Hardy space if 1 (M d ) typically denotes the class of functions with 
the norm 

X>,/iii 

3=0 

where Rj denotes the jth Riesz transform. Here and below we adopt the convention that 
Ro, the 0th Riesz transform, is the identity. This space is invariant under the one-parameter 
family of isotropic dilations, while H l (M. d ) is invariant under dilations of each coordinate sep- 
arately. That is, it is invariant under a t-parameter family of dilations, hence the terminology 
'multiparameter' theory. 

As before, the space i/ 1 (R ci ) has a variety of equivalent norms, in terms of square functions, 
maximal functions and Riesz transforms. For our discussion, the characterization in terms 
of Riesz transforms is most useful: 



ITERATED RIESZ COMMUTATORS: A SIMPLE PROOF OF BOUNDEDNESS 



5 



h 1 (R"*) _ e nn^j*^ 

0<J*<d 



= 1 



R S j s is the Riesz transform computed in the j s ih direction of the sth variable, and the Oth 
Riesz transform is the identity operator. 

4.1. BMO(l^). The dual of the real Hardy space is H X (R S )* = BMO(M*), the t-fold product 
BMO space. It is a Theorem of S.-Y. Chang and R. Fefferman [2] that this space has a 
characterization in terms of a product Carleson measure. 



Define 



(4-2) II&IIhmoo^ = bup l^-'E El< 6 ' 



RCU e*eSig d - 



Here the supremum is taken over all open subsets U cM. d with finite measure, and we use a 



wavelet basis w R . 



4.3. Theorem (Chang-Fefferman BMO). We have the equivalence of norms 

Wf II (if 1 («>*))* ~ II/IIbmo(r <1 )" 
That is, BMO(R«*) is the dual to H 1 ^). 



5. Paraproducts 

We recall a result of Journe [9] , with a recent improvement of Muscalu, Pipher, Tao and 
Thiele [16,17]. (Also see Lacey and Metcalfe [11].) 

Consider the bilinear operators, in fact multi-parameter paraproducts, on functions in M. d , 
given as 

B(/i,/ 2 ) = e*<A>*2></a>*fl>-^b> e R E{-l,+l} 

\K\ 



Rev 3 



we suppress the dependence of this operator on the three indices Ej G {0, l} d , as well as the 
choices of signs e R . 

5.1. Theorem. Recall that d = (d\, . . . ,dt) and that Ej = (£j,i, • • • , £j,t)- If f or all 1 < s < t , 
there is at most one choice of j = 1, 2, 3 with Ej )S = 1, then the operator B satisfies 

B : L p xL q — >L r , Kp, g <oo,± + ± = ±. 

i ft i 'pgr 
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If in addition, E\^\, we will have the estimates 

B : BMO x L p — > L p , 1 < p < oo . 



6. Discrete Analogs of Riesz Transforms 



The Riesz transforms can be obtained from certain types of linear operators which map 
Haar functions to themselves (Haar transforms), a fact demonstrated by Petermichl [18] for 
the Hilbert transform, and Petermichl, Treil and Volberg [19] for Riesz transforms. 

Let us describe the Haar transforms we consider. Fix a dimension d, and a choice of d- 
dimensional dyadic grid T>. Let a : T> — > T> such that 2 d \o~(I)\ = \I\, for all / G T>. Use the 
same notation for a map 



where Sig = {0, l} d — {1} is the set of 'signatures' associated with the wavelets. If <x(e) = 
then := 0. 

Define Q by 



and then extending linearly. There are two facts that we need about these transforms. See 
[19] for the proof of this proposition. 

6.2. Proposition, (a) The operators Q as defined above map L p (M. d ) into itself for all 1 < 
p < oo. 

(b) The Riesz transforms are in the convex hull of the class of operators Q, the convex hull 
taken with respect to the strong operator topology. 

We need tensor products of the operators Q we have just described. For d = (di, d 2 , ■ ■ ■ , d t ), 
and 1 < s < t, let Q s be an operator as above, acting on L 2 (R ds ). We will use the same 
notation for the operator on L 2 (M d ) that equals Q s on the L 2 (M ds ), and is the identity on 
the orthocomplement of L 2 (R. ds ) in L 2 (R. d ). 

6.3. Proposition. The operators 



a : Sig 



SigU{0} 



(6.1) 



def , ff(e) 



(6.4) 



Q = Q! ® ■ • • ® Q t 



extend to bounded linear operators from L p (IR d ) to itself for all 1 < p < oo. 
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Proof of Proposition 6.3. The proof follows by duality. Suppose 1 < p < oo and / G L p (R d )n 
L 2 (R d ) and g G IP' (R d ). Computation gives 

\@f>9)\ = E 



where we have defined the square function 



Hf)(*)= ( E £l</.* 



LR 



Using that 

Haar^ d = {h% : ReVj, e G Sig d -} . 

is a basis for L p (R d ) it is straightforward to show that ||S^(/)||p is an equivalent norm for 
L p (R d ) when 1 < p < oo. This and an application of Cauchy-Schwarz implies 

\(Qf,g)\<\\f\\p\\g\\ P '- 

The above inequality and the density of L 2 (R d ) (~}L p (R d ) in L p (R d ) gives that Q is a bounded 
operator. □ 

7. Proof of Upper Bound for Commutators 

To prove our main Theorem, it suffices to consider the commutators in which the Riesz 
transforms are replaced by choices of the Haar transforms as given in (6.1). 

Define commutators by 



(7.1) C<j(&,/) = [■••[M 6) Q 1 ],-.- ,QJ 

where the Q s are given as in (6.4), and we view the commutator above as acting on L p (R d ). 
In the remainder of this section we prove this proposition. 

7.2. Proposition. The commutators C^ map BMO x LP into LP for 1 < p < oo. 

We make some detailed remarks about the one parameter case. Let Q be as in (6.1), and 
consider the commutator 

(7.3) [M h , , Q}hj = h s X§ - Q hi,h\ 
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There is no contribution if / Pi I' 



0. The other cases yield what follows: 



I cr 

1 = 1' 

I' = a(I) 

V CI and V n cr(I) = 
I' C a(I) . 



(7.4) [AV ; , QM = ^i(i)K(i) ± \I\~ l/2 h7 ( J ) 

±i/r 1/2 < £ ; } ) 
±\a(i)\-^:±\i\^h: 




The first triangular sum corresponding to / C J' is therefore trivial. It illustrates the 
essential cancellation contained in the commutator. It remains to consider the diagonal sums 
1 = 1' and a(I) = I' as well as the other triangular sum /' C a(I). All these do not require 
any cancellation of the commutator as we shall see. 

First, we calculate the diagonal part for 1 = 1' explicitly and obtain the two sums: 



The first sum is essentially a finite combination of operators of the form B(b, Qf) and hence 
bounded. Here we disregard the sign changes as the paraproducts are unconditionally con- 
vergent. We also note the fact that we may change the signatures of the Haar functions, 
which is included in our definition of paraproduct. In the same sense, the second sum has 
terms of the form Q(B(b, /)). The term h\ stems from the product of two Haar functions 
based on I. The signature e may be equal to 1, in which case we still have a convergent 
paraproduct. 

For the other diagonal sum, which corresponds to I' = <j(I), the explicit calculation and 
the reasoning is similar. 

At last, we demonstrate the triangular sum by explicit calculation. We start by summing 
the first summand over /: 



The inner part is equal to YIi-idp X^iW/; ^f)M which on I' equals a renormalized average 

of Qf, namely (Qf, /i}/)|/'| _1//2 . We conclude that the first sum is a combination of para- 
products of the form B(b,Qf). Similarly, the second sum has terms of the form QB(b,f) 
with an indicator falling on /. 



^ ±{b,hi){f,h])\I\' l/2 Q{h £ I ) and Q £ E (6, hf)(J, h})\I\~ 1/2 ^ 




This proves that our commutator is a finite linear combination of terms of the form 

QB(6,/), B(6,Q/) 
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for appropriate choices of Q and paraproducts B. It is essential to note that these para- 
products are of the form that we can permit b e BMO, so that the one parameter form of 
Theorem 5.1 applies to prove Proposition 7.2 in this case. This completes the discussion in 
the one parameter case. 

In passing to the multi-parameter paraproducts, we adopt the notation of Proposition 6.3. 
Thus, we have a choice of d = (d\, . . . ,dt) dimensional dyadic rectangles T>j, and a choice of 
operators Q s , for 1 < s < t, which is an operator as in (6.1) acting on M. ds . We can then 
make an explicit computation of the commutator as in (7.3), namely 

[•••[M, g ,Q 1 ],---Q,] 

The result is a tensor product of the results on the right hand side of (7.4). In other words, 
the multi parameter commutator splits into its components. 

It follows that we can write the commutator as a finite linear combination of terms 

QB(b,f), B(6,Q/) 

for different choices of multi-parameter paraproduct B and different choices of operator Q. 
Thus Proposition 7.2 follows from Theorem 5.1. 
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